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ABSTRACT 


The research described in this thesis deals with an experimental 
method for the measurement of the dual input describing function. -Applica- 
tions of the measured DIDF to actual nonlinear feedback systems will also 
be given. 

The output of the nonlinearity subject to a dual input will contain 
components with frequencies identical to the input and additional harmonics. 
By convention the low frequency input will be called the fundamental signal 
and the high frequency signal the stabilizing signal since the high fre- 
quency signal is usually used to regulate the fundamental output of the 
nonlinearity. 

Two types of stabilizing signals have been used, the sinusoidal 
signal at 100 HZ and the stochastic signal with band width of 55 HZ to 
550 HZ. Each of these signals has been passed through nonlinearities which 
ideally are frequency independent and time invariant. The output was 
measured and converted to the DIDF and then graphed. 

The fundamental output of the nonlinearity subjected to a dual 
input may be measured by summing the output of the nonlinearity with the 
fundamental input of proper sign. Thus if the sum of these signals is 
measured with an RMS meter a minimum may be obtained by varying the funda- 
mental. This eidontin corresponds to the cancellation of the fundamental 
output from the nonlinearity output. If a band pass filter is used around 
the fundamental frequency a sharp minimum may be attained and a DIDF within 


+ 5% of the correct value calculated. 
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i INTRODUCTION 


This thesis is primarily concerned with the measurement of the 
dual input describing functions for both single and multiple valued 
nonlinearities. The method of measurement is applicable to any non- 
linearity but in the interest of simplicity only time invariant and 
frequency independent nonlinearities with three or fewer parameters in 
their input output configurations will be considered. 

To obtain an appreciation of the importance of the DIDF it 
is instructive to outline the development of the describing function. 

The describing function was introduced to approximate the 
action of a nonlinearity in some nonlinear system. Since nonlinear 
elements such as relays, relays with dead band, limiters, limiters with 
dead band, etc., for example, are sturdy, reliable and relatively 
inexpensive, they are widely used in control applications as inexpensive 
amplifiers. As nonlinear elements became incorporated in nonlinear 
systems a method had to be developed to predict the stability of the 
systems. Goldfarb, Tustin, Oppelt and Kochenburger independently 
developed the describing function to provide an approximate analysis 
for the nonlinear systems (Ref. 1). 

The describing function is merely the coefficient of the 
first term of the Fourier series describing the output of the non- 


linearity divided by the amplitude of the fundamental input. The 
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NOMENCLATURE 


transfer function of the low pass filter following 
nonlinearity 


gain constant of the low pass filter 
input to the nonlinearity 
output of the nonlinearity 


amplitude of the fundamental component of the Fourier 


series representing the output of the nonlinearity 


the dual input describing function 

the amplitude of the fundamental sinusoid input to the 
nonlinearity 

the amplitude of the stabilizing signal input to the 
nonlinearity 


the 


standard deviation or RMS value of the stochastic 


input signal to the nonlinearity 


the 


confluent hypergeometric series 


complete elliptical integral of the second kind 


complete elliptical integral of the first kind 


the 


the 


partial differential with respect to A 


ordinary describing function 


parameters used to describe the various nonlinearities 
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describing function is exact only in a system (Fig. 1.1) where an ideal 
low pass filter follows the nonlinearity and thus eliminates all 
fourier components other than the fundamental. Ideally then the 
describing function is the transfer function of the nonlinearity. In 
fact, the describing function ae a limited accuracy. 

Since some nonlinear systems exist which cannot be easily 
controlled by changing the gain without having undesirable side effects, 
it may be desirable to alter the fundamental output by introducing 
another stabilizing signal. Then again, there are some systems which 
already have a secondary signal such as an harmonic signal as in a 
practical feedback system or noise which may be introduced at the input. 
Thus having two inputs a dual input describing function can be defined 
as the ratio of the fundamental output of the nonlinearity to the funda- 
mental input. 

A dual input describing function can not only be used for two 
inputs but it can be extended to the ordinary describing function since 
in ordinary nonlinear feedback systems there are two or more inputs which 
may be significant depending on the low pass filter used. 

The DIDF concept has been used by various authors, among them 
West, Douce and Livesley (Ref. 2) and later Oldenburger and Liu (Ref. 3) 
who introduced an equivalent gain concept very similar to the DIDF. It 
was defined as "the limit of the ratio of the average value of the output 
to the average value of the input as the average value of the input went 
to zero.'' This definition was elaborated on by Boyer (Ref. 4) who intro- 
duced a pseudo describing function which is identical with the dual input 


describing function described by Sridhar (Ref. 5). It was called a 
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FIGURE 1.1 The Exact Describing Function 
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LINEARITY 


w, is the radian frequency of the 
second harmonic of the output of 
the nonlinearity 


w, is the radian frequency of the 
fundamental output of the 
nonlinearity 


6 is an increment in radian 
frequency less than Wo - Wy 
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pseudo describing function because the analytical method used in 

obtaining the function was an approximation to the exact function obtained 
by Sridhar. As it turned out the approximation was very close (identical 
for practical purposes). 

Gibson (Ref. 6) describes a new DIDF which is identical by 
definition to the originally defined DIDF but he extends its use to forced 
systems disturbed by inputs within the operating range of the systems. 
This DIDF is used throughout this thesis. 

This thesis will use a convenient experimental method for the 
determination of the DIDF of the relay, relay with dead band, limiter, 
limiter with dead band, and relay with hysteresis for both sinusoidal and 
stochastic stabilizing signals. In addition the stability of a simple 
nonlinear system (Fig. 1.2) containing the relay with hysteresis will be 
investigated. The nonlinear system will be forced with the two types of 
signals previously mentioned--the sinusoid and the stochastic signal. 
Analytical expressions for the DIDF for both the sinusoid and stochastic 
stabilizing signal exist for the relay, relay with dead band, limiter, and 


limiter with dead band (Ref. 5). 
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EE EXPERIMENTAL METHOD 


a) Representation of the Nonlinearity 


The nonlinearities, other than the relay with hysteresis, 
considered here are time invariant and frequency independent. The relay 
with hysteresis has a slight frequency dependence and the results given 
are not intended to simulate the ideal relay exactly. 

In general the output of the nonlinearity can be analyzed and 
shown to have a fundamental output with the same frequency as the input 
and an infinite number of harmonics with decreasing amplitude. Thus the 


nonlinear device can be described as an input output device as follows: 


NONLINEARITY 
ate “4 oy = £(x) fo y(t) x(t) = INPUT 
aaa 5 y(t) = OUTPUT 


The output y(t) of the nonlinearity may be separated arbi- 


trarily into two components. 


Y(t) = P(t) + H(t) ed) 


P(t) is the fundamental component of the output and has the same fre- 
quency as x(t). H(t) are all higher frequency harmonic outputs of the 


nonlinearity. 
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For the nonlinearity considered here P(t) will at the most be 
a function of five variables--the amplitudes of the inputs plus the three 
parameters of the nonlinearities. In all of the cases two of the parameters 
may fortuitously be eliminated by normalizing the input amplitude with 
respect to parameter. Thus the three remaining variables may be graphed by 


using recurrent plots on a two dimensional graph. 


b) Dual Input Describing Functions for 
Sinusoidal Stabilizing Signals 

The DIDF considered here is the conventional describing function 
adapted for a dual input as first described by J. E. Gibson and R. Sridhar 
(Ref. 6) and incorporated in "Nonlinear Automatic Control" (Ref. 1, page 
420). The DIDF is necessary to obtain some analytical approximation to the 
actual behaviour of a nonlinearity with a dual input. A single sine wave 
of fundamental frequency identical to the fundamental input frequency is 
usually the most logical choice to represent the output of a nonlinearity 
since sine waves can easily be handled mathematically. If this single sine 
wave is chosen as best in a minimum RMS-error sense, then its amplitude is 
the first Fourier coefficient of the fundamental output of the nonlinearity 
element (Ref. 1, page 380). If a dual input exists to the nonlinearity 


then the DIDF may be defined as 
DIDF(A,B) = EAB) (2,2) 


which holds only for a single valued nonlinearity. The dual input to the 


nonlinearity is x(t). 
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= 


x(t) = X,(A) + X,(B) (2.3) 
Xe = A sin at Xe = fundamental signal 
X, ==B sin Dt X, = stabilizing signal 


P(A,B) is the fundamental 
output of the N. L. 


In order to measure the fundamental output we may use the 


following scheme for single valued nonlinearities. 
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The RMS meter reads a minimum when the fundamental output of the 
N. L. is cancelled from y(t). The amplitude of P(A,B) may subsequently be 
read from a - Xo. A better experimental method is to have the needle reach 
its original position. Thus the positional error in the meter is cancelled 
and a - Xe may be divided by two. 


Double valued nonlinearities require a signal 90° out of phase 


with xp which may easily be obtained by passing Xp through an integrator on 
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the analogue computer. 


F_A cos at 8 A cos at 
Xo (B) LOW PASS 
IMadlje SO) — RMS meter 
XE (A) Le 
FILTER 


A sin at (@) . 
a A sin at 
4 


The RMS meter reads a minimum when the real part of the funda- 
mental output of the nonlinearity is cancelled from y(t). Similarly it 
reads a minimum when the imaginary part of the fundamental output is 
cancelled from y(T). The real and imaginary part may then be directly 
measured and the complex dual input describing function obtained. AS 
before a better experimental method is to have the meter return to its 
original position which corresponds to replacing the original fundamental 
output with a signal 180° out of phase. This allows a more accurate read- 


ing to be obtained for the DIDF. 


DIDF (A,B) = eerie + i fe BD. (254) 
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c) Dual Input Describing Functions 


for Stochastic Inputs 


The stabilizing input to the nonlinearity may also be a sto- 


chastic input. 


x(t) = Xp(t) at X, (t) (255) 
Xe = A sin at x = fundamental signal 
Xo = stochastic signal Xo = stabilizing input 


The probability that the amplitude S of the stochastic signal 


X, has a given value S,; is given by the Gaussian Probability-Density 


Function h(x) (Ref. .l7*page 51), 


where P(x) = a‘ exp pose] 
27 0 aor 
HW = average or expected value of the amplitude S 
o = standard deviation, a measure of scatter or dispersion 
o? = variance (second central moment) 
u = 0 for the noise used in this thesis. Note that o the 


standard deviation is also the RMS value of the noise. 

Only noise with Gaussian amplitude distributions is used since 
it is easiest to handle analytically and since it occurs widely in nature. 
Besides, any non-Gaussian noise when passed through a band pass filter 
becomes Gaussian. Consequently Gaussian noise is assured at the output 


of the band pass filter. 
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The fundamental output of the nonlinearity subject to a dual 
input with Gaussian stabilizing signal may be treated in an identical 
Manner to that for a sinusoidal stabilizing signal. The output of the 
nonlinearity will now be a function of the standard deviation o of the 


noise and the fundamental input amplitude A. 


DIDF (A,c) = AE CA,E)) et FRO at (2.6) 
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GEL CONSIDERATION OF THE EXPERIMENTAL TECHNIQUE 


a) General 

To obtain an idea of the validity of the results given in Figures 
Soelte are eORtO pSLilS orandgesial5 toi, LGs6 ter dsttinstruct ive! torexamine-(the 
equipment and its synthesis. 

The Pace 231R analogue computer containing linear passive com- 
ponents of + 0.01% accuracy was used. The error due to physical com- 
ponents was not measurable. However, a true RMS reading meter was used 
for the measurement of the fundamental output and the fundamental input. 
This Briiel and Kjoer Type 2409 electronic voltmeter has an attenuator of 
accuracy better than 1% of full scale deflection. At the 10 HZ where all 
of the output measurements were made a measured accuracy of + 5% was 
obtainable even for a stochastic stabilizing signal. The output fluctu- 
ated especially at low outputs when a stochastic stabilizing signal was 
used. Since analytical results exist for the single valued nonlineari- 
ties considered a quick comparison shows that a maximum error of + 5% 
is easily attained. The close correspondence between analytical results 
and measured results is in itself a check on the ideality of the nonline- 
arity simulations which are considered shortly. 

The errors introduced by the function generating equipment are 
minimal, however they will be mentioned here for the sake of completeness. 
A Hewlett Packard 200CD Wide Range Oscillator was used to generate the 


100 HZ signal. It had a distortion less than 0.5%. The noise was 
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generated by a G. R. random noise generator TYPE 1390-B and was passed 
through a band pass filter (Fig. 3.3) (55-550 HZ) which cut off 60 db per 
decade on the low frequency side. This sharp cutoff was necessary to 
prevent excessive fluctuations in the output of the nonlinearity. 

The measurement technique is shown in Fig. 3.4. All measure- 
ments were done with one meter which was switched back and forth in the 
arrangement shown. Pots Pp and P. were varied to obtain the desired input. 
Switch 4 was used to introduce an imaginary component to the output and to 
the measuring device. P, was used to obtain a minimum reading on the RMS 


4 


Mere. 


b). The Relay 


Simulation of the relay is given in Fig. 3.5. The feedback and 
series limiters sufficiently hardened the output so that there was a net 
change of less than 10 millivolts in the output for a 100 volt change in 
the input. The nonlinearity however was slightly frequency dependent. The 
relay had a finite rise time and settling time which varied as the input 
varied. In general the frequency dependence was not noticeable for inputs 
above 10 millivolts and frequency of 100 HZ or less. 

The error was largest for low inputs, of course, but never 
exceeded the + 5% range. The graphs of the DIDF vs. the amplitude input 
are self explanatory, however they do differ. The slope of the DIDF vs. 
input amplitude graph for a stochastic stabilizing signal is always negative 
while the slope of the DIDF vs. input amplitude graph for the sinusoid 


stabilizing signal has both positive and negative slopes. This phenomenon 
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and its application will be investigated further in a later part of this 


thesis. 


c) The Relay with Dead Band 


The relay with dead band simulation is shown in Fig. 3.5. The 
paths of (+) and (-) incoming signals greater than some predetermined 
minimum, set by the first set of reference pots, is shown. Two amplifiers 
are ganged because the series limiter must end in an impedance of 1 M 
which can only be realized with a gain of 1. The gain of 10 is necessary 
to drive the diodes into saturation completely. The second set of pots 
can be used individually to compensate for voltage drops in the diodes or 
maladjustments of the feedback and series limiters. Since this adjustment 
is rather time consuming it is convenient to have these pots there. 

Upon testing the ideal set up against calculated values it was 
found that the relay did not respond exactly at an amplitude input equal 
to the dead band. At this point the relay output is very sensitive to 
changes in amplitudes of the input and it was found convenient to end the 
dead band a little earlier so that the finite rise time of the relay could 
be partially offset. The relay was again checked and found to correspond 
quite closely to the calculated values. This adjustment is most conveni- 
ently done experimentally by triggering the relays at a lower value of the 
dead band. This correction brings the output of the relay for lower inputs 
closer to calculated values while it worsens the output for higher inputs. 
However, the error is now more evenly distributed over the range of the 


output of the relay. 


. | an 
. “oe 
eld Yo t2sq tetel a ot tedtxyut betagtieeval sd [ftw mokiantiqqe e3t bop 


et 


. ,»eteeds 


or om: a ‘baad basd ditw yeled sat (2 
sit .2,€ .g2% mt owode ef ootsaluate ‘bosd beeb dtlw yslet sfT 


bentmystebs1q smoe cedj teisexg elangte gatmosat (-) bos (+) to efaag 


thee 


‘ezelitiqns owl .mvode at ,atoq sonszsiex io 498 deat? edd yd Isa mumbote 
omit is ceeibteul ns at bos teom rostmil esiise sii seusced begaay 278 
ean at OL to aisg ofT at to oteg es dtiw besiiasx ed yico aso Aoide = 
ajoqg ic “a beuas edit z,yfeislqmos noliayuise otat esbotb sit svigh of 
to asborb ony at aqorb sgaifoy x03 sisansgqities of yi fsubivibnat besev ad aeo 
( Snemseutbs ets ante -arstiatt esiise bus Joadbs9i ed3 to sitoomteolbalem 
| sie ties oiaas svsd ot tnatuovaes al 2: gatewenos sats 19d982 Bk _ 
esau IF seulav beisluotss Jenkngs qu 138 Isebi ef gnkse93 nod 
tines Jugal aierateua SB: led vijsoexs — toa bib etl sit Jad2 bauer 
o3 eviilense ea at jniie ysis1 sdi tatoq ebia tA ‘basd, bask ai7 oF 
edz bite’ o3 jnstneviios bavot esw at bns tuqnt sia Yo esbuttiqas at esgasda ~ 
bLaes gaisaz eel to emis det | sini? ois geds.0a aS ee altiti s beosd besb ‘ 
_baoges:rx0s od ‘baigo baa ‘pateds ntsgs anv tae sit _sa9atte yilstiisg 2 
~tnevaos izom Sk suemsauibe ati “-paplen bedilua iso edt o2 yiseols situp 
a2 io sulav- sewet: a 33 eysisa oct gnizéygtty vd ilesnsabtogze snqb ylins — i 
equgitt tewol ‘203 eee sd3 %0 Jugsuo aris egatrd aoltoarr02 etAT .baed bash’ f ; 
“sesugat aaitgte 308 vugano ada ansatow 2 alice soulav bedaluolss ot aeeala: 


- Yo sanst ods yJavo besudixseib ‘Yaoove 20m wor: “Bk joss sas . 


-.19 - 


FIGURE 3.5 


Analogue Simulation of the Relay 
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Analogue Simulation of the Relay with Deadband 
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A quick check with the scope shows a rise time which varies 
inversely with the amplitude of the input between 5 and 1 ms. A weighted 
mean would be around 2 ms. In order to correct the output for the 2 ms 
delay in triggering one would have to change the dead band of the relay 
as one changed the amplitude of the input. Since this procedure is imprac- 
tical the experimental approach is justified. 

The DIDF for the relay was plotted for both a sinusoidal signal 
and stochastic signal. The sinusoidal plot is in general more accurate 
than the stochastic plot since the stochastic signal had higher frequencies 
which, as noted, the relay cannot respond to. Additionally, of course, the 
meter readings for the stochastic signal fluctuated slightly about an 
average because of the variation in the power output of the noise generator. 


On the whole the readings were within + 5% of the actual calculated values. 


d) The Limiter 

The ideal limiter was reproduced on the analogue computer as 
shown in Fig. 3.14. The set up is quite straightforward and the errors are 
minimal because of its simplicity. 

The results for the limiter were extremely accurate and the 
largest error was due to the measurement technique and the meter used. 

The plots for the sinusoidal and stochastic inputs were quite 


straightforward and need no comment here. 


e) The Limiter with Dead Band 


The simulation of the relay with dead band is shown in Fig. 3.14. 
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The first set of pots set the dead band while the second gave the limiter 
the correct slope and the third attenuated the signal to the proper limit- 
ing value. Again two amplifiers were used in series plus a gain of 10 on 
the amplifier with the feedback limiters to give the necessary gain to 
saturate the diodes even for small inputs. A quick check on the linear part 
of the limiter can be accomplished by bypassing the second pot and every- 
thing following, and feeding into the last amplifier. The output will be 
zero for operating conditions in the linear part of the limiter if the 
slope of the limiter is 1. The procedure is applied separately to each 
branch of the analogue simulation. This procedure shows there is a 
variable error in the output of + 0.7 maximum volts when the limiter is 
operating linearly. This error occurs because the diodes are not com- 
pletely hard and their output varies slightly as they are saturated more 
and more. Consequently for very small outputs of the nonlinearity there 
could be a high percentage error. The results of the measurement tech- 
nique were checked with calculated values. There were no pronounced devi- 
ations although the errors tended to be largest for small values of output. 
On the whole the measured values were within + 5% of the actual values. 

The accuracy of the lower measurements was increased because 
the output tended to be more linear the larger the stabilizing signal used. 
Consequently readings could be taken at larger values and extrapolated to 


the smaller ones. 


f) The Relay with Hysteresis 


The analogue simulation is shown in Fig. 3.19. The four single 
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FIGURE 3.14 Analogue Representation of the Limiter 
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valued nonlinearities were easily set up. The relay with hysteresis how- 
ever could not be made to function properly with an error below + 5%. The 
delay involved in the finite rise times and finite settling times could not 
be adjusted to give a correct output at 10 HZ. Consequently the relay with 
hysteresis was set up ideally for a D.C. input and was treated as a unique 
nonlinearity at 10 HZ. 

After all, the method for measurement described previously is 
applicable to all nonlinearities and a good test of its effectiveness is to 
work with a nonideal relay. 

Using the previously described measuring technique it was found 
that the imaginary component of the output when there was a dual input was 
too small to measure. Consequently only graphs for the real part of the 
DIDF are shown. To check on the validity of the measurements the nonideal 
relay was placed in a feedback loop as in Fig. 1.2 and its stability examined 


experimentally and from the DIDF graphs previously obtained. 
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FIGURE 3.19 


age. 


Analogue Representation of the Relay with Hysteresis 
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IV SOME STABILITY CONSIDERATIONS FOR THE 


RELAY WITH HYSTERESIS G = 5 


In order to obtain an appreciation of the steps involved in 
applying the DIDF and simultaneously checking the applicability of the 
measurements for relay with hysteresis it is instructive to consider the 
following example. 

Consider the relay with hysteresis with G = 5 in a simple 
feedback loop as shown in Fig. 1.2. Proceeding with an analysis in an 
analogous manner as for a single input describing function yields the 
following relationship describing the conditions under which the system 


will oscillate (Ref. 5, page 88). 


; wy athe hen et Kee at 
G (iw) = G (S) = s(S i 10) 2 


DIDF fegt 


Now since the imaginary component of the DIDF was too small to measure 
one may obtain the frequency of oscillation by setting the imaginary part 


of G (iw) equal to zero and solving for w. 


K K 


GE ine) ee 
iw(iw + 10)2 -20w2 + iw(100 - w7) 


7 -20K _ __-i Kw(100 - w*) (4.2) 


400w2 + (100 - w2)2 400wt + w2(100 - w2)2 


388 


~i Kw(100 - w*) 
400wt + w (100 - w2)2 


w= 108rad/sec for £ = 18.59 HZ 


-K -1 
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Additionally |G(iw) | 


K (4.4) 


Thus for a given G(iw) with a known K one can find the DIDF (Eqn. 4.4). 
Knowing the DIDF it is a simple step to go to the appropriate graph which 
represents the DIDF for the particular stabilizing input being used and 
read 2. The procedure is identical whether the stabilizing signal is 
sinusoidal or stochastic. 

Some pertinent results have been tabulated in Table 1 for both 
sinusoidal and stochastic stabilizing signals. The value of K was chosen, 
the DIDF calculated from above and the amplitude of the auto-oscillation 
obtained from the DIDF graphs for various values of stabilizing signal. An 
attempt was made to obtain the value of K at which the auto-oscillation 
would extinguish itself. This value is of course easily attainable from 
the DIDF graphs and occurs either when A is zero, or the value of A at which 
the DIDF peaks and K is obtainable from equation 4.4. The amplitude of A 
was then checked experimentally on the analogue computer. For some DIDF's 
the exact extinguishing point was not attainable and so a range of possible 
values was given. Note the last extinguishing point for a stabilizing 


oC = 6.0 was particularly difficult to obtain. No result was 
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quoted because the DIDF is almost flat at low 4 and hence = can vary broadly 


without affecting the DIDF appreciably. The correspondence of the results 
was better than expected. 

The graph for the DIDF for a dual sinusoid input into a relay with 
hysteresis exhibits definite peaks. These peaks have been investigated 
experimentally and it was found that if for a given DIDF two fundamental 
output amplitudes exist then the lower amplitude is unstable and will 
increase however slowly until the higher amplitude is reached. No attempt 
was made to measure the lower amplitude corresponding to a higher amplitude. 

The instability of amplitudes where the DIDF has a positive slope 
is easily understood if one considers a small perturbation about a point on 
the graph. Thus if the graph has a positive slope and one moves upward one 
increment along the curve, then there is both an increase in the DIDF and 
also the amplitude. The net effect is regenerative and consequently the 
operating point will tend to shift. Similarly, if one moves downward one 
increment then both DIDF and amplitude output are decreased, consequently 
the net effect is to move away from the original operating point. A simi- 
lar analysis of the DIDF curve with negative slope shows that a perturbation 
will tend to cancel itself since either the DIDF is increased and the ampli- 


tude decreased or vice versa. Thus ultimately a stable situation may exist. 
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V AN APPROXIMATE ANALYSIS OF A 


SPECIAL NONLINEAR SYSTEM 


a) The Ideal Relay with Stochastic 

Both variable and fixed gain nonlinear systems exist in industry. 
Thus the output of the former may be controlled by varying the gain while 
the latter is controllable only by introducing an established signal to the 
nonlinearity. If the gain is varied the output cannot be easily predicted 
without using Sridhar's results which require the use of a digital computer, 
or the results which are given in this thesis. In any case there is no 
simple relationship between the two variables. If however the gain K of 
the system is adjusted so that the fundamental output of the nonlinearity 
is quite small and hence usually linear with respect to the fundamental 
input, depending on the value of the stabilizing signal then an approxima- 
tion may be used. 

As an example consider the ideal relay in the nonlinear system of 
Fig. 1.2 stabilized with a Gaussian noise input of standard deviation o. 


The output fundamental amplitude will be given by Ref. 5. 


Lae 
2— Sarg; 2; “Ao (5.1) 
aie 2 20% 


P(A,o) = 


where Fi is the confluent hypergeometric series given by the expression 
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If * is chosen such that 2 £ 2 then if all terms beyond the second in 
the above series are neglected an error not larger than 6.25% will result 


in the DIDF/C. 
Supposing a controllable noise signal existed with a minimum 


standard deviation then the maximum value of A may be obtained which 


“a8 


will make the approximation ae 2 valid. Then a range of DIDF's exists 


(e) 


over which the approximation is valid. 


2 zs DIDF Z Z 
J2n' o . pam o 


Thus a range of values of K for which permissible values of 


(553) 


lw 


DIDF/C exist in the system may be obtained. Hence a K may be chosen 
for a particular system which will allow a simple analytical representation 
of the DIDF in terms of the two inputs. 

In summary then, if there exists a nonlinear system such as in 
Fig. 1.2 a DIDF may be calculated simply from knowledge of the gain K, 
if K is a permissible value as established above. Hence the amplitude 
A of the system output may be calculated simply from the additional 


knowledge of o the standard deviation of the stabilizing signal. 
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Example: 


It is necessary to design a nonlinear system such as in Fig. 1.2 
with an ideal relay which will have an output amplitude controllable by 
a stochastic noise signal which has a minimum value 01, RMS volts. The 


frequency of oscillation of the output may be calculated as on page 37. 


From equation 4.4 it follows that 


DIDF = 


2000 
= (5.4) 


Hence A maximum of the output is a2 Y2 or A maximum is 
A=, {2 peak volts. 


The range of values of permissible gain K may be obtained by 


substituting DIDF = 2000/K for this particular system in the ron 


range 
giving 
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Consequently A may be calculated within + 6.25% in terms of o to give 


A= 2/20 fi. - ts 2000 (5.6) 
1 


This value for A may seem to give a o at which A will extinguish. 
This erroneous impression is caused by the approximation. Actually the 


extinguishing value of o is given by 
Cee pe ol ; 
ae 2 Te for a given value of K (5e4) 


The above expression is easily obtained by setting A = 0 in the original 


expression for DIDF/C given by Eqn. 5.2. 


b) The Ideal Relay with Sinusoidal 
Stabilizing Signal 

The DIDF curves for the ideal relay with a sinusoidal stabilizing 
signal exhibit definite peaks. Operating amplitudes for a simple feedback 
system as in Fig. 1.2 smaller than the amplitude at which the DIDF peaks 
are not attainable because the network is unstable. The DIDF curve has a 
positive slope which causes the amplitude of the output to increase until 
the slope turns negative. 

If the feedback of the network is changed to positive, to attempt 
to compensate for the positive slope, the system output increases or de- 
creases without bound. The operating points on the DIDF vs. A/C curve 
which has a positive slope are not realizable with positive or negative 


feedback. 
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Nevertheless, to check experimental results it is interesting to 
calculate the DIDF when A/C = 0. The normalized DIDF with positive slope 


is given by Ref. 5 as: 


2 
ee eae (Spee 25) 55} BEA (5.8) 
2 1 A ; Ba 


where E &) is the complete elliptical integral of the second kind 
K (2) is the complete elliptical integral of the first kind. 


Set A=0 and by l'Hospital's rule differentiate the numerator and the 
denominator twice with regard to A to obtain the limiting value of DIDF/C 
when A_ goes to zero. 


The limiting value becomes 


DIDF _ _2 


7 an Gey) 


This value of normalized DIDF when A/C = 0 agress for practical purposes 
with the measured value. 

If the system has a negative feedback then self oscillation will 
extinguish at the peak value of the DIDF curve. 

A curve representing the peak values of the DIDF vs. the ampli- 
tude output for the ideal relay with sinusoidal stabilizing signal may be 
calculated as follows. The part of the curve at which the DIDF peaks is 


given by the following equation: 


PIDE - 8 - E(B/A) BSA (Ref. 5) (5.10) 
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Differentiate the equation with respect to A and set it equal to zero to 


solve for the amplitude at which the DIDF/C has a peak value. 
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0 = E(B/A) + a ae (oe1) 


This equation was solved on the digital computer available at the University 


of Alberta to give 


B/A = 0.9089 


The peak values of the DIDF for a given stabilizing signal then occur at a 
fixed amplitude ratio B/A satisfying equation 5.11. The equation of the 


curve passing through the peaks of the DIDF vs. A/C curve now becomes 
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VI AN APPLICATION OF THE DIDF TO THE 


ORDINARY DESCRIBING FUNCTION 


The use of the ordinary describing function for a nonlinear 
autonomous system depends upon the ideality of the low pass filter. To 
check the ideality of any low pass filter and hence the correctness of 
the describing function the DIDF may be used. The second harmonic which 
passes the low pass filter may be used as the stabilizing signal in the 
DIDF. The simplest low pass filter to examine would of course be the 
integrator G (S) = 1/S. However, oscillation can only occur for a double 
valued nonlinearity where the equation G (S) = - ay (where Keq is the 
common describing function) can be satisfied. Since, in practice, the 
imaginary component output of a nonlinearity is usually quite small and 
variable over a small range then the frequency of oscillation of the 
autonomous system would vary. The imaginary component of a nonlinearity 
is usually variable since the practical relay sticks to give it hysteresis. 
The analogue model on the other hand has a finite rise time and settling 
time which affect the imaginary component sufficiently to make calculated 
results impractical. 

A more practical plant G (S) = aa may be chosen in 

Sorrel) 
conjunction with an ideal relay in the configuration of Fig. 1.2 to test 


the common describing function. The frequency of oscillation may easily be 


obtained by setting the imaginary component of G (iw) equal to zero hence 
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Then |G (13) | ate = aa since the lowest harmonic has a fre- 


quency w = 3 rad./sec. Then the amplitude ratio of the two outputs will be 
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where gq is the amplitude ratio of the second harmonic to the first. Now 


the fundamental output of the nonlinearity is given by Ref. 5 as 


P(A,B) = 9 C E(B/A) 


=| 


where A is the fundamental amplitude 
B is the stabilizing amplitude 
P is the fundamental output of the nonlinearity 


E(B/A) is the complete elliptical integral of the second kind 


Now B/A = = and the complete elliptical integral of the second kind may 


be approximated by 


B(B/A) = 3 [1 - 5)? (By) yee 
Then 
eee i eel 
ECA 3) = me C 2 ie 4 (Fe) ° 


Consequently there is less than 0.1% error in the amplitude as predicted 


by the ordinary describing function for this particular example. 
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CONCLUSION 


A simple convenient method for the measurement of dual input 
describing functions was given for any time and frequency independent 
nonlinearity. 

The method was applied to five different nonlinearities; the 
relay, the relay with deadband, the limiter, the limiter with deadband and 
a relay with hysteresis for both a dual sinusoid input and a sine wave plus 
Gaussian noise combination. The results were displayed on graphs. The 
curves of normalized DIDF vs. normalized fundamental amplitude for the 
relay, the relay with deadband, the limiter and the limiter with deadband 
for a dual sine wave input compared favourably with Boyer's results (Ref. 4) 
which were obtained by an approximate analytical method. The curves of 
normalized DIDF vs. normalized fundamental amplitude for the relay, the 
relay with deadband, the limiter, and the limiter with deadband for a 
sine wave input plus a Gaussian noise input compared favourably with the 
analytical expansions obtained by Sridhar (Ref. 5). 

A relay with hysteresis was simulated on the analogue computer 
which approached the ideal kind in accuracy. However, the error for low 
input amplitudes was greater than 10% and consequently the relay was termed 
unique and the DIDF vs. the input amplitude graphed for both a dual sinusoid 
input and a sinusoid plus Gaussian noise input. These unique results were 


tested by comparing experimental results for amplitude and frequency of 


au (ame 


self oscillation in a feedback loop as in Fig. 1.2 to the results obtained 
from the graphs. The comparison as given in Table 1 showed no evident 
discrepancies. 

The imaginary component of the DIDF for the relay with hysteresis 
was not measurable agreeing with the results stated by Boyer that the phase 
shift of the relay with hysteresis is zero (Ref. 4). 

While measuring the amplitude of the fundamental output of the 
nonlinearities mentioned it was noted that for the stochastic plus sinusoid 
input the fundamental output varied very nearly linearly with respect to 
the fundamental input over a given range which could be extended by 
increasing the stochastic signal. This result suggested the approximate 
analysis which was undertaken for the ideal relay. The same type of 
analysis may be undertaken for the relay with deadband, the limiter and 
the limiter with deadband for which analytical results exist. 

In addition the DIDF of the relay, relay with deadband, limiter 
and limiter with deadband can be monotonic with respect to the fundamental 
input for a dual input of sine wave plus Gaussian noise. If then the 
nonlinearity is used in a feedback loop such as in Fig. 1.2 and the 
Gaussian noise is used as the stabilizing signal, then the auto oscillation 
amplitude can be controlled over the complete positive range of values. 

For nonlinearities with a dual sine wave input the auto oscillation cannot 
be similarly controlled. The amplitude may exhibit jump phenomena or 
extinguish at a finite value due to the peaks in the DIDF vs. fundamental 
input amplitude curves. 

The dual input describing function was also applied to an 
example to obtain an appreciation of the accuracy of the ordinary describ- 


ing function. 
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